Photon noise in a random laser amplifier with fluctuating properties 
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We study fluctuations of the number of photocounts measured by an ideal photodetector illu- 
minated by light scattered in an amplifying disordered medium, below the threshold for random 
lasing. We show that the variance of fluctuations and their correlation function carry information 
about fluctuating properties of the medium. A direct link is established between the fluctuations of 
the number of photocounts due to the amplified spontaneous emission (ASE) and the dimensionless 
conductance g of the medium. Our results suggest a possibility of probing amplifying disordered 
media by analyzing statistics of their ASE, without illuminating them from outside by a probe 
beam. 

PACS numbers: 42.55.Zz, 42.50.Ar, 42.50.Lc, 42.50.Nn 



I. INTRODUCTION 



The interplay of multiple scattering and amplifica- 
tion can lead to interesting and still poorly under- 
stood phenomena, culminating in "random lasing" [l|, [3] ■ 
The random laser is a strongly scattering disordered 
medium where losses due to leakage of light through open 
boundaries are compensated by amplification inside the 
medium (sustained due to some external pump mech- 
anism), so that lasing condition is fulfilled despite the 
absence of an external cavity . Similar to conventional 
lasers (i.e. lasers composed of a transparent amplify- 
ing medium inside a high-Q cavity) 'sf, random lasers 
start laser emission only when the pump power reaches a 
threshold. Below the threshold, the system represents a 
"random laser amplifier" : it amplifies light propagating 
inside it (though the amplification is not strong enough 
to compensate for losses) and it is expected to ex- 
hibit properties similar to that of conventional laser am- 
plifiers [6|: increased noise in transmission, amplified 
spontaneous emission, etc. Random laser amplifiers are 
promising systems for studies of Anderson localization of 
light 0, H, @] and, more generally, any interference ef- 
fects in disordered media because, on the one hand, the 
amplification gives increased weight to long-lived quasi- 
modes (or, equivalently, to long wave paths inside the 
medium) and, on the other hand, the system remains lin- 
ear, in contrast to "true" random lasers, operating above 
threshold in strongly nonlinear regimes. However, pre- 
cisely because they operate below laser threshold, laser 
amplifiers require full quantum treatment that properly 
accounts for spontaneous emission processes [1, [a Q- 
One of the quantities for which such a treatment is cru- 
cial is the photocount probability distribution p{n) , giv- 
ing the probability density of counting n photons at the 
output of the amplifier during a given sampling time t. 
In the classical picture, p{n) would be akin to that in 
a medium without amplification, with a larger average 
value of n, but with exactly the same values of normal- 
ized higher statistical moments. In the quantum picture, 
which corresponds to the physical reality observed in ex- 



periments, laser amplifiers are known to add noise to the 
incident light and to generate a noisy light field (am- 
plified spontaneous emission, ASE) even in the absence 
of an external signal [1, |^, Q • For random laser ampli- 
fiers with time-independent properties (quenched disor- 
der) p{n) was studied by Beenakker and Patra in a series 
of papers using the methods of random-matrix 

theory. More recently, Florescu and John have considered 
photon statistics of light emitted by random lasers [l^ . 

In the present paper we study the first two statistical 
moments of p(ji) for light transmitted through or emitted 
by random laser amplifiers with fluctuating properties. 
A simple example of such a system is a suspension of 
scattering centers (particles) in Brownian motion, with 
amplification either in the background medium (typi- 
cally, a laser dye) or in the particles themselves. Such 
media are often used in experiments on random lasers 
[11 [H, m 113, [H, m, Hi . Because of the slow motion of 
scattering centers, all possible configurations of disorder 
are explored continuosly, leading to continuous changes 
of transmission and reflection coefficients of the medium. 
As a consequence, two simple limits exist. First, for 
short sampling times r, well below the correlation time 
tc of, say, transmission coefficients (if the measurement is 
performed in transmission), the measurement yields the 
same result as if the scattering centers were immobile 
(quenched disorder). In the opposite limit of very long 
sampling times t, much in excess of tc, configurational 
average is effectively performed due to the scatterer mo- 
tion and the measurement directly yields the ensemble 
average of the measured quantity (annealed disorder). 
However, only in these two limits a direct link between 
the results obtained for stationary and fiuctuating media 
can be established. In the general case of arbitrary t 
and, in particular, when r ~ tc, the result will depend 
on the particular type of scatterer dynamics and cannot 
be predicted from a priori considerations. It is precisely 
this general situation that we consider in the present pa- 
per. Our analysis is rather general and applies to any 
random laser amplifier with fluctuating properties. To 
be more speciflc, we illustrate our main results on the 
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FIG. 1: A slab-shaped sample of thickness L is filled with an 
amplifying disordered medium represented by discrete, mo- 
bile scattering centers (circles) with amplification of light ei- 
ther in the scatterers, in the surrounding medium, or in both. 
A spatially uniform amplification is achieved through an ex- 
ternal pump mechanism. The slab is illuminated by a probe 
light in the incoming mode a; all other incoming modes a 
and 13 (dashed arrows) are in vacuum states. We are inter- 
ested in the fluctuations of transmitted light in the mode h. 
The polygonal dashed path inside the slab illustrates multiple 
scattering of light with a mean free path £ <^ L. Wavy lines 
with arrows symbolize stimulated and spontaneous emission 
of light by the medium. 



particularly important example of a liquid suspension of 
scattering particles in Brownian motion, where amplifi- 
cation can take place either in the particles, in the liquid, 
or in both. 



II. INPUT-OUTPUT RELATIONS 

We consider a slab of amplifying disordered medium 
with time-dependent properties (see Fig. [T]). The slab 
has thickness L, much exceeding the mean free path £ 
due to scattering, and cross-sectional area A. The mean 
free path is assumed to far exceed the wavelength A of 
light under consideration. At each frequency uj, inci- 
dent and scattered waves can be decomposed over sets 
of 2N{uj) = 4 X k'^A/Air ^ 1 transverse modes sup- 
ported by the slab. Here k = oj/c and c is the speed 
of light in the free space. In counting the number of 
transverse modes we have taken into account 2 propaga- 
tion directions (from left to right and from right to left) 
and 2 polarization states of the electromagnetic field. In 
our notation N{ijj) will stay for the number of left (or 
right) propagating modes at each side of the slab. We 
will use subscripts a and (3 for incoming modes, inci- 
dent on the slab from the left and from the right, re- 



spectively. They will take values a = 1, . . . , N{u>) and 
P = N{ui) + 1, . . . , 2N{lj). The outgoing modes will be 
numbered by 2N{uj) + 1, 3N{uj) for right-propagating 
modes and by 3iV(a;) + 1, ...,47V(w) for left-propagating 
ones. The light incident on the slab is assumed to be 
quasi-monochromatic and corresponds to a single mode 
a = a. It can be in arbitrary quantum state. All other 
incoming modes are assumed to be in vacuum states. 
Our purpose will be the analysis of photon statistics in 
the outgoing mode b on the right of the random medium 
(see Fig. [T]) in the presence of a pump mechanism that 
ensures amplification of light inside the slab. The ampli- 
fication is assumed to be uniform throughout the slab. It 
is characterized by the amplification time ta during which 
the amplitude of a propagating wave is multiplied by e. 
The corresponding amplification length is £a = eta- In 
the context of this paper, we will not need to specify the 
physical mechanism through which pumping is achieved. 

In order to obtain the quantum-mechanical description 
of the problem, we need to quantize the electromagnetic 
field. This can be done in the usual way \^ by replacing 
the amplitudes of the incoming and outgoing modes a(w) 
and their complex conjugates a*(w) by annihilation and 
creation operators d(w) and a^uj), respectively. These 
operators obey the usual bosonic commutation relations, 

[a, {to) , a] [lu')] = 6,jS{u; ~ J) (1) 

where i and j run over all incoming or all outgoing modes. 
Input-output relations between operators corresponding 
to incoming and outgoing modes in the case of an ampli- 
fying disordered slab with time-i ndep endent properties 
have been estabhshed in Refs. [13, [Hi [13, [HI : 

CLb (w) = ^ tab (w) CLa {i^) + ^ T pb (w) CLp (w) 

Q=l I3=N{lj) + 1 

+ Y.v;,{u:)c\{u^) (2) 

7 

Here tab and r^b are transmission and reflection coeffi- 
cients from the incoming modes a and /?, respectively, 
to the outgoing mode h. The last sum in Eq. ^ de- 
scribes spontaneous emission of light in the amplifying 
medium, with 7 running over "objects" (e.g., atoms or 
molecules) emitting at frequency lo. To be precise, the 
phenomenological model leading to Eq. Q is based on 
the assumption that the interaction of light with a lin- 
ear amplifying medium can be modeled by coupling of 
optical modes to a bath of inverted harmonic oscillators 
[23 | . The number of terms in the last sum of Eq. ^ is 
then the number of oscillators in the bath, whereas the 
bosonic operators obey 

[c,{uj),c]{J)]^5,,5{uj-J) (3) 
{cMc]{oj')) = -S,j5{oj - cj'Ulu) (4) 

with ri{LL!) = [exp{haj/kT) — 1]^^ and negative effective 
temperature T < 0. 
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In a medium with time- dependent properties the input- 
output relations ^ should be generalized and become 



Of, (uj) 



dio' 



tab (W, w') CLa [uj') 

a=l 



2N{lj') 

ri3b{uJ,uj')dp{uj') 



(5) 



Modes at different frequencies are not independent any- 
more: the incident field at a given frequency u' can give 
rise to outgoing waves at different frequencies uj lo' . 
The physical meaning of this can be understood on the 
simple example of a suspension of scattering particles. 
Because of the thermal motion of particles, scattered 
light will suffer from Doppler broadening and hence will 
not be monochromatic even if the incident light were. 
The Doppler broadening of transmitted, reflected, and 
spontaneously emitted light are described by the coeffi- 
cients tabi(^,(^'), fpi,{uj,uj') and respectively. 

Let us now show that Eq. (O indeed describes a 
medium with time-dependent properties. To this end 
we consider a realistic situation of a measurement device 
sensitive only to light within a narrow frequency band 
Aw ^ wq around the frequency loq of the incident wave. 
The bandwidth Aoj can be due either to some technical 
limitations of the device or to a band-pass filter put in 
front of it on purpose. Taking spectral filtering into ac- 
count, the result of any measurement performed on the 
outgoing mode h can be expressed through an inverse 
windowed Fourier transform cn,{t) of (Xb{uj)'. 



ab{t) 



1 



uJo+Auj/2 



due 



*ab(w) 



(6) 



ljo~Alj/2 



An expression for db{t) can be obtained from Eq. ([5]) by 
applying the inverse windowed Fourier transform to its 
right-hand side. For the first sum of Eq. ([5]), for example, 
this gives 



a)o+Aa)/2 



due' 



oc 



N{u)') 

dui' y tab {lU,Uj') da {io') 



oo 



(7) 



We now introduce sum and difference variables u) = 
(lu + Lo') /2 and = u — lo' and define a new transmis- 
sion coefficient depending on these variables: tab{<u, fi) = 
tabiUjUj'). As a function of il, the latter will be assumed 
to be different from zero only in a narrow frequency band 
Su; <^ Aw around = 0, which corresponds to slow vari- 
ation of the properties of the random medium on the 
time scale of 1/Aw. As a consequence, not only the in- 
tegration over LU, but also the integration over uj' will be 



effectively limited to the pass band Aw of the filter. We 
can therefore replace N{uj') in the upper limit of summa- 
tion in Eq. ([7]) by A^(wo). To lighten the notation we will 
drop the argument wq of N(ujq) in the following and will 
write N = N{ujo). Next, we assume that as a function 
of LU, tab vary significantly only on the scale of F ^ Aw. 
Then, because w, w' and hence w remain within the pass 
band of the filter, we can replace the first argument w of 
tab by Wo and obtain 



1 N /.a)o+Aa)/2 



dw'iab (wo, w — w') (w') (8) 



The next step is to represent t^b {'-^o, lu — lu') as a Fourier 
transform of tabiwo,t') and to perform the integration 
over LU. This yields 



N 



oo poo 



-iuj't' —iujQ (t — t') 



J2 dw' dt't 

J -oo J-oo 

tab{uJ0,t')da{uj')5^^{t-t') (9) 



where we defined 







/ dwe"* 




'-Auj/2 


Aw 


( Aluu\ 




smc 




V 2 ; 



(10) 



and sinc(x) — sin(a;)/a;. S^^{u) tends to the Dirac delta- 
function in the limit of Aw — > oo. Because tab i^Ojt') is 
a slowly varying function of its second argument (time) 
and only vary significantly on the scale of 1/^w 1/Aw, 
while SaljIu) is peaked around t = t' and has a width 
1/Aw, we can replace t' by t in the argument oitab (uq, t') 
and take the latter out of the integral in Eq. ([9]). The 
remaining integrals over t' and w' are easily calculated, 
leaving us with 



N 



^'^tab (wo,i) da (t) 



(11) 



where da (t) is the windowed inverse Fourier transform 
of da {lu) defined as in Eq. ([6]). Finally, because the three 
remaining sums in Eq. ([5]) can be treated in the same way, 
we obtain the input-output relations in the time-domain: 



N 



2N 



Of) (t) = X! + X! ^1^'' 

a=l P=N+1 



J2v;b{t)c\{t) 



(12) 



where tab {t) = \J1-Ktab (wo,t) and similarly for Vjib and 
w*^, except that the Fourier transforms are defined with 
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FIG. 2: Frequency scales of the problem. The wide dashed 
curve schematizes the slow dependence of the transmission 
coefficient tahiy,'^') on lj = u}' . The narrow solid curve de- 
picts the dependence of this coefficient on w ^ uj' . When 
uo' is varied, the latter curve shifts along uj axis, keeping its 
overall shape and remaining centered at uj' . Its maximum 
touches the curve tctb{i^, The width 5uj of the solid curve is 
much smaller than the width F of the dashed curve. A band- 
pass spectral filter limits any measurement to a frequency 
band centered at ujq. The pass band of the filter Alu obeys 
Slj <^ Auj <t. r. 



opposite signs of i for v^h and 



1 



2tt 

2tT Jujg-Auj/2 



dne-''^*v;,{u;o,n) (13) 
diue-^'^'cUij) (14) 



We will use the input-output relations (jl'ip in the remain- 
der of the paper. They are valid under the assumption 
of Suj ^ Au! ^ r, where Scu is the largest of the widths 
of Fourier transforms of tab (t) , r^f, (t) and u*^ (t) ; Auj 
is the width of the frequency band to which the mea- 
surement is limited, and F is the typical scale of varia- 
tions of tab (wo, t), rpb ((^0, and u*^ {ojq, t) as functions 
of the carrier frequency wq- Physically, the first condition 
Slu <^ Auj means that we are considering the limit of slow 
variations of the properties of disordered medium (small 
Su and hence long typical times ~ l/Suj), such that the 
broadening of the spectrum due to these variations is well 
below Aw. The second condition Aw <C F allows us to 
neglect variations of the frequency-dependent properties 
of the medium (such as the index of refraction, the mean 
free path, etc.) within the pass band Aw. This signifi- 
cantly simplifies calculations and makes the results more 
general and not limited to a particular random medium. 
The different time scales of the problem and their hier- 
archy are illustrated in Fig. [2l 

Working in the time domain requires commutation re- 
lation for time-dependent operators. These are easily 
obtained from Eqs. ([1]) and ([3]): 



[a, (t) , a] it')] = S,,Sa^ it - t') e— o(t-f, (^5) 



-iLO{-){t~t') 

[c, [t) , c\ it')] = <5,,(5a. {t - t') e»-«(*-*') (16) 



Equation Q in time domain becomes 

(q (t,) c\ (t,)) = ~5,,r^{ujo)SAUt - t')e*""^*-*'^ 



(17) 



where we have taken into account the condition Aw <ti 
wq and approximated r]{uj) by ?/(wo) for all w e [wq — 
Aw/2, Wo + Aw/2]. Note that r/(wo) < in an amplifying 
medium. 

It will appear convenient in what follows to define in- 
tensity transmission, reflection, and spontaneous emis- 
sion coefficients: Tab{t) = \tab{t)\'^, R0b{t) = |r-/3b(t)P 
and V^-yb(t) — |u-yb(t)p. Requiring that the commutation 
relations (fT5|) hold for the outgoing mode 6 and mak- 
ing use of the slowness of variations of tab{t), r/^bit) and 
Vybit) on the scale of 1/Aw, we obtain an energy conser- 
vation law for our problem: 



Tb (t) + Rb (t) - Vb (t) - 1 



(18) 



where we defined the total transmission, reflection, and 
spontaneous emission coefficients; Tb{t) — '^a'^°'b{t), 
Mt) = EpRpbit) and Vbit) = EyVibit). We note 
here that these definitions of the total transmission and 
reflection coefficients Tb and Rb are equivalent to those 
usually used in mesoscopic wave physics (23l . [2^ , except 
that they are defined for a given outgoing mode (mode b) 
with summations in their definitions running over all in- 
coming modes, whereas the usual definition is for a given 
incoming mode with summations running over all outgo- 
ing modes. For a reciprocal medium, to which we will 
limit our consideration here, these two definitions are of 
course equivalent. We also note the appearance of a new 
quantity Vb{t) in Eq. (fT5|) . These quantity vanishes in the 
absence of amplification, when Tb (t) + Rb (t) = I. In an 
amplifying medium, its appearance in Eq. (jlSp illustrates 
that amplification [i.e. the condition Tb {t) + Rb (t) > 1] 
necessarily implies additional noise [the last sum in Eq. 
P^ ] . This can be interpreted as a manifestation of the 
fluctuation-dissipation theorem in this particular case. 



III. AVERAGE NUMBER OF PHOTOCOUNTS 

We now consider an ideal photodetector illuminated by 
the outgoing mode b. The operator corresponding to the 
number of photocounts registered by the photodetector 
during a sampling time r is 



fib = dt al{t)ab{t) 



(19) 



Using Eqs. ^T7\ and the fact that all incoming 

modes except a single mode a are in vacuum states, one 
readily calculates the quantum mechanical expectation 
value of hb'. 



{rib) = 



dt Tab {t) {di [t) da {t))-Tj{LJn)^ / dt Vb (t) 

Jo 

(20) 
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FIG. 3: Average amplified spontaneous emission coefficient 
Vh as a function of the ratio of the slab thiclcness L to the 
amplification length La for 3 different thicknesses L/l = 10, 
100 and 1000. Dashed fines show asymptotes for x — > and 
a; — + TT in the case of L/l — 10. 



We will use angular brackets (. . .) to denote quantum 
mechanical expectation values from here on. 

Equation (|20p was obtained for arbitrary but given 
time-dependent functions Tab{i) and Vb{i). (hb) can fluc- 
tuate from one realization of these two random functions 
to another. We can perform ensemble averaging, denoted 
by a horizontal bar (...), over all possible random real- 
izations of Tab{t) and Vb{t): 



{rib 



Tabifla) 



tAuj 
2tt 



ri{LL>o)Vb 



(21) 



where {fia) is the average number of photocounts that 
would be counted by the same photodetector in the in- 
cident mode a during the same time t, and the random 
processes Tab(t) and Vb{t) are assumed stationary, so that 
their statistical moments are independent of time. 

Equation ((2T|l contains two terms. The first one de- 
scribes transmission of the incident light through the slab 
of amplifying medium, whereas the second one repre- 
sents the amplified spontaneous emission (ASE) of the 
random medium. The average transmission and sponta- 
neous emission coefficients Tab and Vb entering Eq. (1^ 
can be calculated in the diffusion approximation (see Ap- 
pendix El : 



Vb - 



1 sin{£/La) 
N sm{L/La) 
sinje/ La) + sin[(L - 
sm{L/La) 



- 1 



(22) 
(23) 



where La = \filaji is the macroscopic amplification 
length and la is the inverse amplification coefficient. 
Both Tab and Vb diverge at LjLa — tt which we iden- 
tify as a laser threshold, following previous literature 
[HlillljIlSjIlH- Obviously, our approach based on the lin- 
ear equations (|12p can only be applied below the thresh- 
old, i.e. for LjLa < tt, which is a condition for having a 



random laser amplifier and not a random laser. The av- 
erage transmission coefficient of a random laser amplifier 
Tab was previously studied in Refs. [13, [3l|. However, 
the spontaneous emission coefficient Vb was not intro- 
duced previously, even though ASE of a random ampli- 
fying medium was studied in Refs. [l^, [13 • We show Vb 
as a function of L/La in Fig. [31 The fraction of photo- 
counts due to amplified spontaneous emission in the total 
signal can be found by dividing the second term in Eq. 
(ED) hylf^: 



-T]{uJo)VbTAoj/{27r) 

(fib) 



(24) 



We will find this quantity useful later. It varies between 
(no spontaneous emission) and 1 (no incoming light). 

We will see that in the following an important role will 
be played by the average dimensionless conductance of 
the disordered slab in the absence of amplification [La — > 
oo), equal to the transmittance T and defined as 



3^ 3 L 



(25) 



Because we do not treat the dependence of Tab on a and 
b in full (in reality. Tab is not independent of a and b but 
contains slowly varying dependences on the directions of 
propagation of corresponding modes [2^), we added a 
factor I in front of the sum over modes in Eq. . This 
makes the final expression for g consistent with its stan- 
dard definition [23, [IJ] ■ A transition from diffuse mul- 
tiple scattering to Anderson localization is expected at 
5 « 1 [sj . In the present paper we only consider the case 
of g > 1. 



IV. VARIANCE OF PHOTOCOUNT 
DISTRIBUTION 

We now proceed to the calculation of the fiuctuations 
in the photocount number. These fluctuations can be of 
classical or quantum origin. Classical fluctuations arise 
because the transmission coefficient Tab{t) fluctuates in 
time and from one realization of disorder to another. The 
fluctuations of Tab(t) can be characterized by its autocor- 
relation function 



CTabT,, (t) 



Tab [t') Tab [f + t) 



r. 



(26) 



ab 



Quantum fiuctuations appear due to the nonclassical na- 
ture of light in the incoming mode a and spontaneous 
emission of light by the atoms in the random medium. 
The degree of "quantumness" of the incident light can be 
quantified by the Mandel parameter Qa = var(na)/ (fia) — 
1 [6], where var(na) = ("a) ~ ('^a)^, whereas the inten- 
sity of spontaneous emission is controlled by the Bose- 
Einstein function ri{ujQ). Because both the incident light 
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and the spontaneous emission have to cross the whole 
or a part of the disordered medium in order to be de- 
tected in the outgoing mode b, quantum fluctuations will 
be "weighted" by the fluctuations of classical coefficients 
Tab and Vf, in the final result for the variance of fib. In 
particular, we will see that, in addition to ([26]) . the fol- 
lowing correlation functions appear: 



Tab {t')Vbit' + t) 

(^T^bVAt) = ^ 



TabVb 



Vb (tp Vb [f + t) 

vl 



(27) 
(28) 



The calculation of correlation functions CT^^Ta^{t)^ 
^TabVbit) and CvtVb{t) is presented in the Appendices 
IB | [Cl and ID | respectively. In order to calculate these cor- 
relation functions, we excluded the immediate vicinity of 
the laser threshold from our consideration and assumed 
that A = 1 — L/irLa » l/y^, where, we remind, g ^ I- 

With these definitions in hand, we start by calculating 
the quantum mechanical expectation value of hi , 



{hl)= dt, dt2{al{h)ab{h)dl{t2)db{t2)) (29) 



In addition to the assumption of slow variations of tab (t) , 
rpbit) and w*f,(0 on the scale of l/Aw, adopted in the 
previous sections, we restrict our consideration to sta- 
tionary, monochromatic incident light in the mode a, for 
which J, = {di{t)da{t)) and Fa = (fit (i)at (t')a,(t)a,(t')) 
are time-independent and {d\{uj)da{i^')) — 2TTlaS{uj — 
uj')S{uj — loq). a relatively simple result for (n^) can 
be then obtained in the limit of long measurement time 
T > l/Aw: 



inl) 



{hb)+Fa / dt / dt'Tab{t)Tab{t') 

Jo Jo 



2Iarj{ujo) 



A f-T f-T 

— dt dt'Tab{t)Vb{t') 

27r 7o Jo 



dtTab{t)Vb{t) 



V i^o) 



Alu 
~^ Jo 



dtVbit) 



^ fdtvm 



(30) 



We will characterize fluctuations of hb by the their vari- 
ance var(nf,) = {hi) — {{fib))'^ , divided by the square of the 
average number of photocounts: 



SI 



va.T{hb)/{{hb)y 
;r an ensemble c 
realizations of the disordered medium, use Eq. (|2ip for 



To calculate , we average Eq. (|30p over an ensemble of 



and definition p4| for ip. This yields 



-5? 



oT. (^) 



{rib) ^ \NAujJ 



2(l + C^T„.v.(0)) 



1-ip 



(l + CvbvM) 

1 + 3 



Qa 



TabTa, 



(31) 



where rp 



= {2It) /^(l - t/T)CT^^T^,{t)dt and sim- 
ilarly for (5|, ^ and Sy^ y^ . This equation is the main 
result of the paper. 

To discuss the meaning of different terms in Eq. (|3ip . 
let us first consider the situation when the spontaneous 
emission is absent {ip = 0). Eq. (ISDl th en reduces to the 
result obtained previously in Ref. [2^]. It contains the 
shot noise l/(f^b) due to the fact that the energy of the 
electromagnetic field in the outgoing mode b can be ab- 
sorbed by the detector only in discrete portions (quanta) 
huj, the classical noise due to the fluctuations 

of the classical transmission coefficient Tab, and a com- 
bined term Qa(l + S^ h)/{^a) that vanishes for light 
in the coherent ("classical") state for which Qa — 0, and 
can be either positive or negative for nonclassical light 
(e.g., light in Fock or squeezed states) [1^. When the 
spontaneous emission is present (tp > 0), the situation is 
more complex and all terms of Eq. (j3ip have to be taken 
into account. The terms of Eq. PT|) containing Cy^Vb(O) 
and Sy^y^ originate from the spontaneous emission alone, 
whereas those containing Ct^^ vi,(0) and S^^^y^ represent 
the interference of the incoming light with the sponta- 
neous emission. 

Because Eq. (I5T]) is quite lengthy and difficult to ana- 
lyze in its general form, we will restrict our consideration 
to two particularly interesting cases. The first case cor- 
responds to a strong incoming wave such that la ^ and 
the fraction of spontaneous emission in the total number 
of photocounts in the outgoing mode b is small: ip ^ 1. 
In this limit, we will expand Eq. (j3ip in power series in 
p and will keep only linear terms (Section IIV A[) . The 
second case corresponds to the absence of incoming wave 
{la — and </?=!). The outgoing light in the mode b is 
then due uniquely to the amplified spontaneous emission 
(Section [IVB|. 



A. Strong incident wave 

When the incident wave in the incoming mode a is 
strong and the fraction of photocounts due to the am- 
plified spontaneous emission in the outgoing mode b is 
small, we can develop Eq. ([3T|) in series in p <^ 1 and 
keep only linear terms. To be specific, we assume that 
the incident wave is in the coherent state {Qa = 0). After 
a substitution of Eq. (j24p for p, the normalized variance 
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FIG. 4: Correlation function Ct„,v,(0) = TahVb / {T abV i,) - 1 
as a function of the ratio of slab thickness to the amplification 
length, X = L/La- Dashed lines show asymptotes for a; ^ 
and X —> TT. 

of photocount fluctuations can be written as 

61 = =i={l-2,7(^o)n[l + CT„,y,(0)]} 
{nb) 

27jiiOo)Vb fNAu;\ . . . . 

+ ^Lt.. (32) 

For the amplifying medium with time- independent, static 
properties, S^^^j,^^ = CT„tT„,,(0) ~ 1 and S^^^y^ = 
Cr„(,Vi,(0) ^ 1/g. The latter quantity is shown in Fig. 
|4] as a function of the ratio of slab thickness L to the 
amplification length La- Equation (|32p then becomes 

<5,2 = =l={l-27j{cjo)Vb[l + CT^,vM]} 
{nb} 

+ Ct„,t„.(0) (33) 

This result differs from the result of Patra and Beenakker 
because averaging over realizations of disorder was 
performed in a different way in that work. A result con- 
sistent with Ref. [l^l is obtained by setting CT^^TabiO) = 

and Ct^.vAO) = «: = [1 - 2r]iuJo)Vb]/^. Obvi- 
ously, this corresponds to neglecting fluctuations of Tab 
and Vb from one realization of disorder to another. 

Surprisingly, the different way of averaging over dis- 
order may have an important influence on seemingly 
very general conclusions that one makes by analyzing 
the variance of fib. Indeed, in the limit of weak am- 
plification I^f, ^ 1 and for g 3> 1 and tAlu ^ 1 , we 
readily obtain from Eq. ([55]) that varn^ — S^Hfib))^ — 
{fib) -\- {{fib)Y + 2T]{LjQ)Vb{fib){TAuj/2n), where the first 
term describes the shot noise, the second — large inten- 
sity fluctuations due to disorder in the medium [26l | , and 



FIG. 5: Normalized variance of photon number fluctuations 
in transmission of a strong incident wave through a slab of 
amplifying disordered med ium. For this figure, we set slab 
thickness L/l = 100, '{hb)^ = 10, 2nIa/NAio = 10, g = 100 
and ri{uJo) = — 1 (i.e. T 0~). The two curves correspond 
to L/La = (no amplification) and 1. Note that the curve 
corresponding to L/La ~ 1 is slightly below the L/La = 
curve. 



the last one can be interpreted as an excess noise due to 
the amplified spontaneous emission. Note that because 
ri{tjjQ) < in an amplifying medium, the excess noise 
is negative. This result, also illustrated in Fig. [51 is a 
consequence of the interplay between quantum fluctua- 
tions, due to the interference of the incident light with 
the amplified spontaneous emission, and classical fluctu- 
ations, due to the randomness of the transmission and 
spontaneous emission coefficients Tab and Vb- It seems 
to be in conflict with the positive excess noise found in 
Ref. fl^, but the solution of this contradiction lies in 
different ways of averaging over disorder: we consider 
varfih = (n^) — {{fib))"^ instead of var'fif, = (n^) — {fib)'^ 
in Ref. When we calculate the latter quantity, we 

find var'rif, = {fib) — 2r](ujo)Vb{ fib) and the excess noise 
is positive, in agreement with [12l |. An important con- 
clusion of this analysis is that if theoretical results are 
compared to experiments, care should be taken to en- 
sure that averaging is performed in an appropriate way. 
In a medium with fluctuating properties, the simultane- 
ous quantum and disorder averages (. . .) that we use is 
the only experimentally feasible option because both are 
typically replaced by time averages and it is hence not 
possible to perform the quantum average (. . .) without 
performing averaging over disorder (. . .) [25l[27]|. 

Let us now analyze the behavior of Eq. ([5^ as a func- 
tion of {fib). To be more specific, we restrict our con- 
sideration to the case when the random medium repre- 
sents itself an ensemble of point-like scattering centers in 
Brownian motion with diffusion coefRcient D^- In this 
situation, the motion of scatterers introduces a charac- 
teristic time tc = ^to{£/L)'^, where to = l/4fc^DB- This 
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1 2 
Amplification x = L/La 

FIG. 6: Variance of fluctuations of tlie spontaneous emission 
coefficient 14: Cvf,Vt{Q) ~ '^b f^b — 1 as a function of tlie 
ratio of slab thickness L to the amplification length La . 



S 




Time y = \/t/tc 

FIG. 7: Correlation function Cv^v^it) as a function of the 
square root of normalized time y — \/TJt'c for several values 
of the ratio L/ La- 



is the typical correlation time of transmitted light. We 
find convenient to introduce the average number of pho- 
tons {fib) ^ detected during a sampling time equal to tc- 
Three distinct regimes can then be identified in the de- 
pendence of 5l on {hbj- First, for (rib) < 1, Eq. ([5^ 
is dominated by the first term on its right-hand side: 
^ ^/{fib}- Second, when (hb) becomes of the order of 
1, another two terms on the right-hand side of Eq. ((32|) 
come into play. In the limit of {fib) > {fib)c^ find from 
the results of Appendix iBl that ^T^^T^t ^ ^i(^6)c/ + 
{C2/g){{fib)c/{f'-b))^^^, where Ci and C2 are constants 
and the two contributions to rp are due to the cor- 

relation functions C'tHt^jCO ^^'^ respectively. 
These are the short- and long-range correlation functions 
defined in Appendix [Bj Using results of Appendix El 
we also find S^^^y^ oc {1/ g){{fib) ^/ {hb))^^"^ ■ As long as 
(fit,) < g'^{fib)^, the first term of i5|. ^ dominates and we 

obtain 6^ oc l/{hb), similarly to the case of (nb) < (fib)c, 
but with a different coefficient. However, the truly large- 
{fib) behavior is reached for (fib) > g'^{'n,b)c^ when the 
second term of S'^^^T^b ^TotVt become larger than the 

first term of (^T^^T^i,' this hmit. We 

illustrate the three regimes in the dependence of b\ on 
{fib) and transitions between them in Fig. O 



B. Noise in the amplified spontaneous emission 

We saw in the previous section that the interference 
of a strong incident wave in the coherent state with the 
weak amplified spontaneous emission in an amplifying 
random medium modifies noise properties of the trans- 
mitted light. We now consider another special case of the 
general Eq. pTjl : a situation when no light is sent into the 
random medium from the outside (i.e. (fia) = 0) and the 



signal measured in the outgoing mode h is due uniquely 
to the amplified spontaneous emission of the medium it- 
self. One readily sees from Eqs. pTjl and ((3T|l that in this 
case 

^ = -^?7(^o)n (34) 

= 7= [1 - ^(^o)n (1 + CvbVb (0))] + (35) 
{nb) 

Once again, we recover the result obtained previously 

[13], SI = i/W) + 2tt/t/:^uj, if we neglect fluctuations 

of Vb from one realization of disorder to another and set 
CvbVbi^) = and Sy^Vb =0- It is true that these fluctua- 
tions are weak: CvbVbi^) ~ 1/5 <C 1 and 5vbVb 1/g 
in the limit of weak ampliflcation, but they grow when 
the laser threshold is approached and even diverge at the 
threshold (see Figs. Eland [5]). 

In order to study the behavior of Eq. ((35l) as a function 
of (fib), we first consider how Sy^^y^ depends on the sam- 
pling time r. 5y^y^ is an integral of CvbVb{^) niultiplied 
by (1 — t/r) and divided by t/2. For an ensemble of scat- 
tering centers in Brownian motion, Cv'^Vb(i) is calculated 
in the Appendix |D1 and we show it in Fig. |7l for several 
values of L/La- In the limit of short sampling times t, 
we have 5'^^y^ = CvbVbi^) with CvbVbi^) given by Eq. 
(|D3p of Appendix [D1 and shown as a function of L/La in 
Fig. m In the limit of long sampling times, the scaling of 
5y y with T is determined by the long-time asymptotics 
oiCvbVbit)- It follows from Eq. (|D2| of Appendix that 
CvbVb{t) = \/tcJt/9 with tc = |io(^/i)^ in the limit of 
t > tc- This yields 5y^y^^ = {^/ig)yJtjT. This asymp- 
totic behavior of 5y^y^ is illustrated in Fig. |8l Because 
according to Eq. ([M)) . (fib) oc t, we can replace the ratio 
tc/r by a ratio {fib) ^1 {fib) with (fib)^ denoting the aver- 
age number of photocounts in a sampling time equal to 

tj,. We therefore obtain 5y^y^ — {8/3g)J {fib) ^1 {f'-b)- 
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FIG. 8: The most interesting contribution to the normalized 
variance of photocount fluctuations, Sy^y^, as a function of 
the square root of normalized sampling time y = y/rjil for 
several values of the ratio L/La. 



We are now in a position to analyze the behavior 
of 51 as a function of {fib)- At small {fih), ^VbVi, '^^-n 
be neglected with respect to the first term in Eq. ([35|) 
and we obtain S'^ cx l/{hb). But even though initially 

1/2 

small, SvbVi decays with (rib) only as I /{fib) , which is 
due to the long-range nature of the correlation function 
C^VbVbit)- As a consequence, Sy^y^ dominates the results 
in the limit of large 



[nb): 




{rib 



(36) 



The transition between the small- and large- (n;,) regimes 
in the dependence of on (hb) takes place when 
l/(nb) ^ '^Vi.Vi,' i-e- at (hb) g if (fib)^ > 5 or at 
{fib) ^ 9^ / {fib) c otherwise. We illustrate the transition 
between the two regimes in Fig. [3] 

Our result ()36|) calls for two important comments. 
First, it makes an interesting link between random lasers 
and Anderson localization 3- Indeed, Eq. states 
that in the limit of large average photon numbers, the 
noise of the amplified spontaneous emission is propor- 
tional to 1/g, with g — the dimensionless conductance 
of the disordered medium — determining the closeness of 
the Anderson localization transition, expected at g — 1 
Q. This suggests that 5 of a disordered sample can be 
measured by observing fluctuations of its amplified spon- 
taneous emission. Another interesting observation is that 
Eq. (|36p contains the average number of photocounts in 
a time interval equal to the correlation time of hght scat- 
tered in the disordered medium, t^. Therefore, it appears 
possible to obtain useful information about the dynamics 
of scattering centers in a random medium (e.g., the diffu- 
sion coefhcient of particles in suspension) from the fluc- 
tuations of amplified spontaneous emission. This "noise 
spectroscopy" can be made even more efficient by using 
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FIG. 9: Normalized variance of photocount fluctuations 5^ 
resulting from the amplified spontaneous emission of a dis- 
ordered slab of thickness L. The three curves correspond 
to three different values of the ratio of L to the a mplifi cation 
length La- Other parameters are set to L/l = 100, {fib)^ ~ 10, 
g = 100 and r]{uo) = —1. 



the full equation (|35p to fit the whole curve S'^{{fib)) and 
not only its large- (rib) part. 



V. PHOTON CORRELATION SPECTROSCOPY 
IN A RANDOM LASER AMPLIFIER 

In addition to the variance, another important quan- 
tity characterizing a fluctuating signal is the autocorrela- 
tion function of its fluctuations. The technique that uses 
the autocorrelation function of photon number fluctua- 
tions to characterize disordered media is known as pho- 
ton correlation spectroscopy (PCS) HI]. Typically, one 
measures the autocorrelation function of photocounts de- 
tected by a photodetector in a mode b of light scattered 
by a disordered sample: 



where 



C„(,„(,(t) 



fibit') = 



{fib{t')hb{f + 1)) 



{fib{t')) ■ {fibit' + t)) 



1 



t'+T 



dt"al {t") ab in 



(37) 



(38) 



is the operator corresponding to the number of photo- 
counts measured by the photodetector in the time inter- 
val [t',t' + t]. In a medium with fluctuating properties 
(like, e.g., a suspension of dielectric particles), Cmmif) 
contains information about the type and the intensity of 
fluctuations intrinsic to the medium (like, e.g., the type of 
particle motion — Brownian, directed, arrested dynam- 
ics, etc., — and its typical velocity). In optically dense 
media where light undergoes multiple scattering, PCS 
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is most often called difFusing-wave spectroscopy (DWS) 

MM- 

The autocorrelation function ((37|) of light transmitted 
through a slab of amplifying random medium can be cal- 
culated using methods developed in Section HVl In fact, 
the normalized variance of photon number fluctuations 6^ 
is nothing else than C„^„^(0). In the present section, we 
will restrict our consideration to t > t, i.e. we will be in- 
terested in the correlation of photocounts corresponding 
to two nonoverlapping time intervals of equal duration 
T. Using Eq. commutation relations ([15]), pB]) and 
Eq. (fT7|) . we readily obtain 

t' + T t'+t + T 

{ni,{t')nb{t' + t)) = J dh J dt2 
t' t'+t 



Tab{h)Tab{h){ai{ti)a\{t2)aa{ti)aa{t2)) 



2^7j{uo)Tabiti)Vb{t2){ai{h)aa{h)) 

ZTT 



+ ( ^ ) v'MVbih)Vbit2) 



(39) 



where the limit of long sampling times r ^ l/Aui 
was taken. Note that here, in contrast to Eq. of 
Section [Wl we did not assume that {al{ti)aa{ti)) and 
{al^{ti)a'l{t2)aa{ti)aa{t2)) are time- independent. We will 
not need this assumption in this section. 

We now average Eq. over realizations of disor- 
der and assume that the correlation functions CT^jTai, (i), 
Cx^^Vtit) and Cy^Vbii) do not vary significantly on the 
time scale of r. (Please note that this situation is dif- 
ferent from the one considered in Section IIVI where the 
main idea was precisely to take into account fluctuations 
of Tab and Vb on times shorter than r.) After a series of 
algebraic transformations we finally obtain 



Cmrif, it) 



(l-^)2{[l + C„„„„(t)]CT„,T„.(t) 

C„„„Ji)} + 2^(l-^)CT„,n(t) 
^^Cv,v,{t) (40) 



where is the fraction of photocounts due to the ampli- 
fied spontaneous emission defined by Eq. the cor- 
relation function of photocounts in the incident mode a, 
C„^„^(<), is defined in the same way as C„^„^(t), and 
CT„i,T„i,(i), Cr^tVbit) and Cy,y,(t) are defined by Eqs. 
(HlD-dMl). 

The structure of Eq. (gO]) is similar to that of Eq. (PT|) . 
The first term is the only to survive in the absence of 
amplified spontaneous emission {(p = 0). In this case, 
Eq. (PP]) reduces to the result previously obtained in Ref. 
|2a |. The second term, proportional to Lp{l — ip), repre- 
sents the interference of external light with the amplified 
spontaneous emission and requires that both are present 
(i.e. that < (p < 1). Expanding Eq. in series in 
ip 1, we can see that, similarly, to what happens to 6^, 



the amplified spontaneous emission provides a small cor- 
rection to ip — result. The most interesting limit of Eq. 
dJO]) is the limit of ip = 1, corresponding to the absence 
of external illumination. In this limit, only the last term 
of Eq. (|40)) . representing the autocorrelation function of 
the amplified spontaneous emission, survives and 



(41) 



For an ensemble of point-like scattering centers in Brow- 
nian motion, this correlation function is given by Eq. 
(ID2p of Appendix ID] and we show it in Fig. [71 For weak 
amplification x — L/La ^ 1, it is small in magnitude 
1/g) for 5 3> 1, but it grows with x. Cv^Vbi^) di- 
verges at the laser threshold x = it (see Fig. [S|). Inter- 
estingly, Cvi,Vh{t) is long-range in time and decays only 
as 1/ Vt. Long-range spatial correlations of similar origin 
have been previously reported by Patra and Beenakker 
[lll | . A new feature of our result ([4T|) is that it contains 
information about the dynamics of fluctuations in the 
disordered medium (e.g., information about dynamics of 
scattering centers). Measurements of Cn^n^{t) can there- 
fore be used for spectroscopy of disordered media, even 
in the absence of external light source. 



VI. CONCLUSIONS 

We presented a theoretical study of noise in the number 
of photocounts measured by an ideal, fast photodetector 
illuminated by a single mode b of light emerging from 
an amplifying disordered medium (random laser ampli- 
fier) with fluctuating properties. Assuming that the pho- 
todetector is sensitive only to light in a small frequency 
band Aw, we derived general expressions for the normal- 
ized variance and autocorrelation function Cn^n^it) 
of photocount fluctuations in the limit of long sampling 
times r ^ l/Auj. Because light in the mode b contains 
a mixture of transmitted incident light and light sponta- 
neously emitted and then amplified by the medium, the 
fraction of light due to the amplified spontaneous emis- 
sion ip appears as a natural parameter that controls the 
behavior of 5^ and C„j„j(t). For Sf, two distinct regimes 
were identified and studied in detail. First, when <IC 1, 
the overall behavior of remains qualitatively similar to 
that a,t ip — 0. For a suspension of scattering centers in 
Brownian motion, for example, we found cx l/{hi,) in 

1/2 

the limit of (hi,) and S"^ oc \/g{fn,) in the limit 
of (rib) —t oo, with g the dimensionless conductance of 
the disordered sample in the absence of amplification. 
Second, when 95 = 1, i.e. in the absence of external il- 
lumination, our result describes statistical properties of 

light spontaneously emitted by a random laser amplifier. 

1/2 

We have found 5^ cx l/g{nb) in the limit of (fib) — > 00 
for this case. 

For Cn^mit), the two regimes <C 1 and cp = 1 exist as 
well, with the most interesting result Cn^nti^) oc \/g\/t 
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obtained for the latter one. This long-range time correla- 
tion of the amplified spontaneous emission could possibly 
be used to probe dynamics of scattering centers in ran- 
dom media. 

Our results allow us to establish an interesting link 
between random amplifying media and Anderson local- 
ization. We have found that for the amplified sponta- 
neous emission, both and C„j„^(t) are inversely pro- 
portional to the dimensionless conductance g of the dis- 
ordered medium (we remind that the localization tran- 
sition is expected at g = 1). Our results suggest that 
measurements of and Cn^nt{t) in ^ random laser am- 
plifier could allow a precise determination of g. 
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and the solution of Eq. (jAip reads 



sinh 




sinh 


^q^ + xHL-z^) 


\J + sinh 







(A2) 



with r — {p,z), z< = min(z,z') and z> = max(z,z'). 
The average intensity of light inside a slab of amplifying 
medium /^(r), resulting from a plane wave incident on 
the slab in the mode a, obeys the diffusion equation (|Aip 
with x^ = and —S{z — i)/A on the right-hand 

side. Here La = \fUaJi is the macroscopic amplification 
length and Xjla is the amplification coefhcient. /^(r) can 
be therefore expressed through the Green's function (|A2p 
as 



This work is supported by the French ANR (project 
No. 06-BLAN-0096 CAROL) and the French Ministry of 
Education and Research (Research-Educational Network 
"Quantum optics of random media"). 



/„(r) ~ -- / dVG(r,r')(5(z'-£) (A3) 



with the integral over the volume of the slab. The average 
transmission and reflection coefhcients from the incoming 
modes a and /3, respectively, to the outgoing mode h are 



APPENDIX A: AVERAGE TRANSMISSION, 
REFLECTION AND SPONTANEOUS EMISSION 
COEFFICIENTS OF AN AMPLIFYING 
DISORDERED MEDIUM 



^i:Dd- 1 sin(£/Lj 

Tnh = — — /a(r) Iz^L = -77 . . (A4) 



dz 

AttD d ^ 



NsiXi{L/La) 



k^ di 



1 sin[(L-^)/L,] 
Ip{v)U^,=- ^.^^^^^^^ (A5) 



In this Appendix, as well as in the three next Appen- 
dicies [Bl [C] and [Dl we will ignore the vector nature of 
electromagnetic waves, except for the expression of the 
number of transverse modes TV — fc^A/27r which is cal- 
culated with a proper account for two independent po- 
larizations states. Such a simplification is sufficient for 
our purposes because the key properties of average trans- 
mission, reflection and spontaneous emission coefficients, 
as well as the properties of their correlations functions, 
that we use in the main text, are the same for vector and 
scalar waves. 

For distances exceeding the mean free path £, the trans- 
port of light in disordered media can be described in 
the diffusion approximation. A fundamental quantity 
through which many other quantities can be expressed 
is then the Green's function of the diffusion equation 
G(r,r') obeying 



D{W'-X') G(r,r') 



6{r- 



(Al) 



where D — c£/3 is the photon diffusion coefficient, c is 
the speed of light in the medium, and x is a parameter 
that we will specify later. In a slab of thickness L and 
surface A ^ L^, perpendicular to the z axis, boundary 
conditions are G = at the surfaces z = and z = L, 



Within the approximations that we made in this paper, 
these results are independent of the mode indices a, (3 
and b. More accurate considerations show that both T„ h 
and Rfjb exhibit slow dependences on a, (3 and h [23l.[24|. 
In addition, a coherent backscattering cone exists in re- 
flection [11]. These corrections to T^b and R^b might 
be important if a quantitative comparison with experi- 
ments is attempted, but we neglect them here because 
they introduce additional technical complications in the 
analysis without bringing any qualitatively new physics. 

The total transmission and reflections coefflcients, as 
well as the spontaneous emission coefficient then follow 
from Eqs. (|^ and ((J5|) : 



Tb = Tab — — 



sin(£/La) 



sin(L/La) 



- ^- sin[(£-^)/£„] 
^ fim{L/La) 

Vb = Tb + Rb-l 

sm{l/La)+sm[{L-t)lLa] 



(A6) 
(A7) 

1 (A8) 



sin(L/La) 

This completes derivation of Eqs. (|^^ and (|23p of the 
main text. 
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APPENDIX B: CORRELATION FUNCTION OF 



The correlation function CT^^T^^{t) defined in Eq. ([26]) 
can be decomposed into a sum of a large but short- 
range part cl^^'^j, ^ (t) and a small but long-range part 

^T^bT^.it) SHlail,!!!]. Because we will need it later, 
we consider a more general object CT^^T^,^{t) where a' 
can be different from a. The short-range part of the cor- 



relation function, C 



(1) 



(i), can be obtained from the 



autocorrelation function Caa' (r, t) = {Ea{r, t')E^,{r, t' + 
t)), where Ea{r,t) is the complex electric field inside the 
disordered medium illuminated by a wave in the incom- 
ing mode a. If the time dependence of Ea{r, t) is due to 
the motion of scattering centers in the medium, Caa' (r, i) 
obeys the same diffusion equation as the average inten- 
sity /Q(r) but with = +7^(0 ~ l/-^aj where 
= ~ Qq' is the difference of transverse compo- 
nents q of wavevectors k = (q, k^) of the incoming modes 
a and a', ^{t) describes decorrelation due to the motion 
of scattering centers and depends on the type and inten- 
sity of motion. For Brownian motion with a diffusion 
coefBcient DB,j^it)^ 3t/2toi^ with to = I/AU'^Db- Us- 
ing the circular Gaussian statistics of Ea (r, t) , we obtain 



^ ab ^ ' 



d 



Caa' (r, t) 1 2 



3inh(V72(t) - 1/Ll) sm{L/La) 



sm 



(Bl) 



This correlation function is short-range in both time and 
space. It decays exponentially with time t and it vanishes 
for a ^ a' . In the absence of amplification {La ^ oo), 
the correlation time of the fluctuations of Tab is tc — 
\to{llLf. 

In order to calculate the long-range part of the cor- 



(2) 



(i), we apply the well-known 



relation function, C 
Langevin approach which we generalize to include both 
amplification (as in Ref. i30i. i31i]) and motion of scatter- 
ing centers (as in Ref. [33|). The fluctuation of intensity 
(5/Q(r,t) — Ia{v,t) — /a(r) obeys a diffusion equation 

D{W^ + l/Ll)6Ia{r,t) = divjjr,0 (B2) 

with uncorrelated external Langevin currents 



jL'\r, t')£}(r, t' + t)^ S.,^ \Caa' (r, S{r ~ r') 



27r^c2 



|2 r/ r\ 



(B3) 



Using the relation 



we obtain 



^T„,(iO^Wt' + t) = ^(^) A 
d d r'^ d ~ d ~ 



dz" 



X \Caa'{v"Ml=.'=l 



(B5) 



where 6Tab{t) = Tab(t) — Tab and G(z,z") = 
/ (Pp"G{r, r") is the q = Fourier transform of G(r, r") 
with respect to p — p" . Evaluating the integral in Eq. 

2 

(jB5|) . dividing it by T^f^, and taking the limit of L/£^ 1, 
we arrive at the following result: 



^tX,W = -F2{L/La,L^nt) + Acil) (B6) 



where g = ^N£/L is the dimensionless condutance and 



Axy^ sinh^ \/ — x'^ 



sinh 2 \/y^~~x^ 



2xy'^ ~ {y"^ ~ a;^) sin2x 



(B7) 



For y — 0, corresponding to t = 0, Eq. (|B7[) becomes 

cotan X 1 



F2{x,y) 



sin^ X 



(B8) 



It is quite generally known that another contribution to 
the long-range correlation function C'^^^j. ^ ^ [t) exists (23l . 
[23 |. By symmetry, this contribution depends on Aq^ = 

— qf,, in the same way as Eq. (|B6p depends on Aq^. 
In our case, b — b' and Aq^ — 0. Therefore, the final 
expression for the long-range correlation is 



C 



(2) 

To.bT, 



F2{L/Lan{t)L) 



+ F2{L/La,WjHt) + Aql) (B9) 



The correlation function (|B9p is long-range in space (it 
is of the same order for a — a' and a a') and time (it 
decays only as l/Vt in the long-time limit). Because of 
the prefactor I/5 <C 1, it is, however, small in magnitude, 
at least as long as the amplification is weak (a; = L/ La <^ 
!)■ 



APPENDIX C: CORRELATION FUNCTION OF 

T^b{t) AND Vb{t) 



STabit) = J dPp^JIa(v,t)U=L (B4) 



The correlation function of Tab{t) and Vi,{t) de- 
fined by Eq. (|27p requires calculation of an average 
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STab{t')SVt{t' + t). Because of Eq. ^ and definitions 
of Ti,{t) and Rb{t), this average can be represented as 



Cll). we obtain 



ST^bit')6Vbit' + t) 

= 6T^bit')ST^'bit' + t)+J2 ST^bit')SRfjbit' + t) 

a' p 

= TabYTa'bCT^^T^,^(t) + Tqfc RfSbCr^^Rg^ {t) 
a' /3 

(CI) 

Let us first consider the first sum. In the absence of am- 
phfication (i.e. for La — > oo), C'Tq(,t^/[, (0 is known to 
have three distinct contributions denoted as C'^^\ C^^^ 
and C(3) ill. The first two of these were ex- 

phcitly calculated in Appendix [Bl C'^' is of order 1, 
but it vanishes for a' ^ a. C^^' is a factor 1/g smaller, 
but it is of similar mae; nitude for all a'. Finally, C^^) 
is yet another factor 1/g smaller and, similarly to C*^^-', 
does not vanish for a' ^ a. We assume that this hierar- 
chy of correlation functions does not change significantly 
in the presence of weak amplification (i.e. for La finite 
but still La ^ L), which is indeed the case for C*^^^ and 
as can be seen from the results of Appendix |B] 
Then, when the different types of correlation functions 
are substituted into the first sum of Eq. (IClll and the 



2 2 

summation is performed, they will yield T^^, {N/g)T^fj 

and {N/g'^)T^ij, respectively. Because in this paper we 
consider a situation when N/g ^ L/i ^ 1 and g ^ I, 
the contribution of C*^^^ will dominate the result and the 
two other contributions can be neglected. 

Suppose now that amplification is strong and we are 
close to the laser threshold (i.e. that A = 1 — L/irLa ^ 
1). The behavior of equal-time correlation functions C^'-* 
in this situation has been studied by Burkov and Zyuzin 
[H: C(i) ~ S^^,, ~ l/gA2 and ~ l/g^A'^. 

Now C*^^-* will dominate the first sum of Eq. (|C1[) only if 
1/ yAg, otherwise the largest contribution will come 
from C^^-* . In the present paper we assume that the con- 
dition A ^ ^/y/g is fulfilled, which also ensures that 
C'^' is always much smaller than C*-^-*. We can there- 
fore neglect all contributions to the correlation function 
CT^tT^,,{t) in Eq. other than C^^^ for both weak 

and strong amplification, excluding only a narrow region 
in the vicinity of the laser threshold which requires a sep- 
arate treatment and will not be considered here. Among 
the two terms contributing to cl^^^j. [see Eq. (jB9P ]. the 
first one F2{L / La,^{t)L) will give the main contribution 
to the sum over a' in Eq. (jCip because it is independent 
of Aq^, while the second term decays with Ag^. We will 
neglect the second term of Eq. (jB9p in the summation of 
Eq. (lUTj) . 

Obviously, a similar reasoning applies to the second 
sum of Eq. (jCip . except that C^^^ correlation is always 
zero. CT^f,Rf3i{t) can be obtained along the same lines 
as CTabTndh Appendix [b1 Inserting the first term 

of Eq. (|B9j) and a similar result for Ct^^r^^^^ it) into Eq. 



CT^,v,{t) - - X |-^cotan- 
coth \J ip- — ( 

sin X 



2 sinh^ \/y'^ — x"^ 



1 - 



xy^ 



(C2) 



where x — L/La and y = ^{t)L. In particular, for t — Q 
Eq. ([02]) becomes 



g 2x sin x 
3 

— cos X ^ 1 



sinx V 2 



cos X 



(C3) 



APPENDIX D: AUTOCORRELATION 
FUNCTION OF Vb{t) 

Similarly to Appendix [Cl we represent 
6Vb{t')5Vb{t' + t) as 

5Vb{t')5Vb{t' + t) = J2STc.bit')ST^'bit' + t) 

a, a' 

+ Y.^^Pbit')SRp>bit' + t) 

13,13' 



+ 2j2ST^b{t')5Rpb{t' + t) (Dl) 

For the same reasons and under the same conditions as 
in Appendix [cl the main contributions to the three sums 
of Eq. (jPip come from C^^^ correlation functions. As 
compared to Eq. (jCip . Eq. (|Dip contains a new corre- 
lation function SRfjb{t')5Rp'b(t' + t) which, however, can 
be calculated in the same way as the two others. When 
divided by vl, Eq. (|DT|1 yields 



Cnvt,(0 



1 



1 



g 4a;?;2y^j2rr^sin2 (I) 
X |2a; [j/^ -|- (x^ — y"^) cosx] cotanh \/ y^ ~ x"^ 



sinh^ V "'^ — 



2 I J2 



2xy' 



[x? " ly^ + x^ cosh(2 \J ip- ~ x^ ) 1 sm 2: 



} (D2) 



For f = wc obtain 



1 



Cv.vM = - X 



9 16a; (sir 



[4x(2 -I- cos x) 



1 sin a; — 4 sin 2x + sin 3.t1 



(D3) 
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